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Abstract 

^ There are recent considerations to increase the photomultipher density in the IceCube de- 
Q-i tector array beyond that of DeepCore, which will lead to a lower detector energy threshold 
pH and a relatively huge fiducial mass for the neutrino detection. This initiative is known as 
"Precision IceCube Next Generation Upgrade" (PINGU). We discuss the possibility to send 
psj a neutrino beam from one of the major accelerator laboratories in the Northern hemisphere 
^ to such a detector. Such an experiment would be unique in the sense that it would be the 
only neutrino beam where the baseline crosses the Earth's outer core. We study the de- 
0\ tector requirements for a beta beam, a neutrino factory beam, and a superbeam, where we 
consider the cases of both small ^13 and large 6*13, as suggested by the recent T2K and Dou- 
2 ble Chooz results. We illustrate that a fiavor-clean beta beam best suits the requirements 
^ of such a detector, in particular, that PINGU may replace a magic baseline detector for 
small values of 613 - even in the absence of any energy resolution capability. For large 613, 
, ^ however, a single-baseline beta beam experiment cannot compete if it is constrained by the 
X CERN-SPS. For a neutrino factory, because of the missing charge identification possibility 
in the detector, a very good energy resolution is required. If this can be achieved, especially 
a low energy neutrino factory, which does not suffer from the tau contamination, may be an 
interesting option for large 61^. For the superbeam, where we consider the LBNE beam as a 
reference, electron neutrino flavor identiflcation and statistics are two of the primary limita- 
tions. Finally, we demonstrate that in principle the neutrino factory and superbeam options 
may measure the density of the Earth's core at a sub percent level for sin^ 26*13 ^ 0.01. 



^Eniail : j t ang@physik . uni-wuerzburg. de 
^ Email : winter @physik. uni-wuerzburg. de 



1 Introduction 



Neutrino oscillation experiments have provided compelling evidence that the weakly inter- 
acting neutrinos are massive particles [l]. Given three generations of massive neutrinos, 
there must be two characteristic mass squared splittings (Am|^ , Am2i) and three mix- 
ing angles {Ou , 6'i3 , 6*23) as well as a Dirac-type CP violation phase (5cp) affecting neu- 
trino oscillations. Disappearance of muon neutrinos, which is mainly driven by |Am|^| 
and 6*23, has been observed in atmospheric neutrino oscillation experiments, such as Super- 
Kamiokande [2], and in the MINOS long baseline experiment [sj. Disappearance of electron 
neutrinos has been observed from solar neutrino oscillation experiments relating to 612 [4], 
whereas Am2i has been strongly constrained by the KamLAND long baseline reactor neu- 
trino experiment [sj. The CHOOZ short-baseline reactor neutrino oscillation experiment [6] 
has provided a limit sin^ 2^13 < 0.1. Recently there has been a 2.5a hint for non-zero ^13 
from the T2K experiment [T] and the Double Chooz experiment [s] , see Refs. |9l|l0| for global 
fits. Since it is too early to claim a ^13 discovery, we leave 61^ as unknown parameter in this 
study. We quahtatively distinguish the "large 6*13" (6'i3 discovered by the next generation 
of experiments) and "small ^13" {613 constrained by the next generation of experiments) 



cases. These cases are roughly separated by sin^ 2^13 = 0.01 [ll,12 . The most important 
questions are then: is Amgj^ > (normal ordering) or Amg^^ < (inverted ordering), what is 
the value of ^13, and is there leptonic CP violation (CPV). Other open questions concern the 
possible existence of light sterile neutrinos, as indicated by recent MiniBOONE antineutrino 
results 13 confirming the LSND anomaly 14 , and a (possibly energy- or environment-) de- 



pendent superluminal motion of neutrinos 15 . Many of these questions can most likely not 
be addressed with existing experiments 12 , including the reactor experiments searching for 
^^13. A next generation of experiments, using accelerator-based technology in the Northern 
hemisphere, has been studied in international efforts, including the International Neutrino 



Factory and Superbeam Scoping Study 16 -18 , the ongoing International Design Study 



for the Neutrino Factory (IDS-NF) UM, and the Euronu design study "A High Intensity 



Neutrino Oscillation Facility in Europe" 20 



In the Southern hemisphere, the IceCube Neutrino Observatory 21 has recently completed 
deployment at the South Pole, including an array with denser string configuration called 
DeepCore with a lower threshold (9(10) GeV. Given that this technology is very well 
established and the cost for future potential upgrades is relatively well predictable, there 
is recent activity for the extension of IceCube in the low energy direction, called PINGU 
(Precision IceCube Next Generation Upgrade). The primary physics applications for this 
experiment are those similar to that of other megaton-size water Cherenkov (WC) detec- 
tors, such as proton decay, detection of supernova neutrinos, and detection of atmospheric 
neutrinos. However, such a detector could also conceivably receive a neutrino beam, as dis- 
cussed for some superbeam-upgrade and beta beam proposals in the Northern hemisphere 
(for example, CERN-Frejus or J- PARC to Hyper-Kamiokande) . For the first step of such an 
upgrade, 20 additional strings are discussed leading to a fiducial mass of 1 Mt at a few GeV, 
increasing to about 10 Mt at 10 GeV. In a future phase, for which the detector parameters 
are not yet specified, there may be even a substantially larger fiducial mass. Although the 
baseline from any of the major accelerator laboratories in the Northern hemisphere is known 
to be suboptimal for neutrino oscillation studies, in particular for CP violation, we study in 
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this work if the large increase of statistics by this relatively huge detector mass may compen- 
sate for the loss of statistics from the long baseline. More specifically, consider a beta beam, 
neutrino factory beam, and a superbeam, and establish the detector requirements for each 
of these beam classes for small and large 613. Sec. |4] discusses our detector parameterization. 
These requirements may serve as guidance for the optimization of the detector, if it were 
intended to receive a neutrino beam. Our primary focus is the optimization with respect to 
the most often used performance indicators, including the ^13, mass hierarchy (MH), and 
CP violation (CPV) discovery reaches as a function of (true) 6*13 and 6cp- Note that the 
neutrino oscillation physics using such a long core-crossing basehne is quite unique because 



of the "parametric enhancement" 22 , 23 of the oscillation probability; for a detailed dis- 
cussion, see Sec. |2} Conversely, the density of the Earth's core can be probed with such a 
neutrino beam [24], which we study in Sec. [7j Earlier works studying a neutrino beam to 



the South Pole include Refs. 25 26 



As far as the different beam classes are concerned, we distinguish setups for small 613, where 
only an upper bound on 6*13 can be expected from the next generation of superbeam and 
reactor experiments, and large ^13, where ^13 will be discovered by the next generation of 
experiments. The setups for small ^13 should be able to access as small as possible 6*13, 
whereas the setups for large ^13 should have an optimal Sqp reach. For a neutrino factory, 
the small and large ^13 setups are specified in the IDS-NF Interim Design Report [27]. For 
small values of di^, a two-baseline setup with Li ^ 3000-5000 km (for CP violation) and the 



so-called "magic baseline" 28 with L2 — 7500 km (for degeneracy resolution and the mass 
hierarchy measurement) is proposed, with a muon energy i?^ = 25 GeV. For large values of 
^13, a single-baseline setup with a smaller muon energy and shorter baselines is preferable, 
where the exact parameters depend on the detector properties. For a magnetized iron 



detector, ~ 10 GeV and L ~ 2200 km has been obtained in Ref. 29 . For a beta beam, 



qualitatively similar setups have been established for small 30 and large 31 values of ^13, 
using two and one baselines, respectively. For a superbeam, typically only single-baseline 
options are considered in the literature since the intrinsic contamination of the beam with 
electron neutrinos limits the ^13 reach. An example of this is the Long-Baseline Neutrino 
Experiment (LBNE) wide band beam experiment from the Fermi National Accelerator 



Laboratory (FNAL, Batavia, Illinois) to Homestake (North Dakota) in the US 32 . In this 
study, we consider modifications of the above beam setups. For small values of ^13, we test 
if the magic baseline detector for a beta beam or neutrino factory could be replaced by 
PINGU; see Sec. |5} For large values of 6^13, we consider single-baseline options using a beta 
beam, neutrino factory beam, or superbeam; see Sec. [6j Details on the beams, detector 
requirements, and reference setups can be found in Sec. |3j 



2 Oscillation physics with a core-crossing baseline 

The relevant information on sin^ 2^13, 6cp, and the mass hierarchy at long-baseline neutrino 
oscillation experiments is typically inferred from the appearance probability Pg^t = Pve^y 
(or P^e)- For a constant density profile, it is typically appropriate to expand this probability 
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to second order in sin 2^13 and the hierarchy parameter a = Amgi/Amg^^ ~ 0.03 as 33 35p 



P ~ 
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Here A31 = Aml^L/{4:E), A = 2V2GFneE/Aml^ with rig the electron density in Earth 
matter, and the sign of the second term refers to neutrinos (positive) or anti-neutrinos 
(negative). Note that the sign of A also depends on considering neutrinos or anti-neutrinos. 
The T- inverted probability, P^e, may be obtained from Eq. ([T]) by changing only the sign of 
the second term. 

Especially if sin26'i3 ~ a, i.e., all the four terms in Eq. ([T]) contribute equally, the correla- 
tions and degeneracies affect the extraction of the individual parameters 36 . A convenient 
method is to combine two baselines 27,30 : a shorter baseline sensitive to all four terms in 
Eq. (|T|) t o measure CP violation (in the second and third terms), and a longer "magic" base- 



line 



28 



7500km for degeneracy resolution. At this second baseline sin(y4A3i) = 0, 
meaning that only the first term in Eq. survives allowing for a clean sin^ 2^13 and 
mass hierarchy measurement. The condition sin(ylA3i) = can be fulfilled independent of 
oscillation parameters and neutrino energy, and depends only on the matter density pro- 
file. However, if two baselines are combined, a wider baseline range produces a similar 
effect 
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The large fiducial mass of PINGU may compensate for any disadvantage 
from a longer baseline, such as lower event rates and some re-occurring dependence on 6cp 
(which may prove useful for CPV studies). Therefore, we will study two-baseline setups 
where the second baseline is replaced by that to PINGU. For large values of ^13, the first 
term in Eq. ([T]) dominates, and is the primary background for the CP violation measure- 
ment. Therefore, systematic uncertainties are important, but their detailed investigation 
goes beyond this study. Note, however, that broader beam spectra p59^], or longer base- 
lines that include more oscillation peaks 40 tend to be less sensitive to systematics than 



narrow-band off-axis beams operated at the first oscillation maximum. 

Another peculiarity obtained from Eq. ([l]) is the effect of the MSW resonance A — )■ 1. 
In this case the first term will be minimally affected by the l/L^ decrease of the flux for 
a wide baseline range 



41 



meaning that large event rates can be expected even for very 
long baselines. This condition, however, is energy-dependent, and is satisfied for the MSW 
resonance energy 



E ~ 



13200 cos 2^ 



13 



p [g cm" 



(2) 



"'^ Although this expansion is not particularly accurate for very long baselines 35 , especially for baselines 
crossing the Earth's core, it is sufficient to illustrate our main qualitative points. 
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Figure 1: Left panel: Matter density profile for a baseline L = 11810 km for the Reference Earth 
Model [ 42| (thin curve) and a three-step approximation (thick curve) . Right panel: P^,^ appearance proba- 
bility for these two matter density profiles (same line styles; normal hierarchy). 



for neutrinos and the normal mass ordering, or antineutrinos and the inverted mass ordering. 
For the average density of the Earth's mantle and small sin^2^i3, one finds i^rcs — 7GeV, 
while the average density of the Earth's (outer) core gives i^res — 3 GeV. 

Since the major high energy accelerator laboratories are in the northern hemisphere, some 
potential baselines to PINGU will cross the Earth's core. More specifically, some baselines 
ar43 



• CERN-PINGU: 11810 km 

• FNAL-PINGU: 11620 km 

• JPARC-PINGU: 11370 km 

• RAL-PINGU: 12020 km 



Considering oscillation physics, these baselines are very similar and the physics outcome 
will be the same. Therefore GERN-PINGU: 11810 km is chosen as an example for this 
study. We show in Fig. [T| left panel, the matter density profile for L = 11810 km for the 
Reference Earth Model (thin curve), and a three-step approximation (thick curve) consisting 
of the average densities in the qualitatively different layers. Compared to Eq. ([T]), there are 
two very different matter densities present at the same time, leading to a "parametric 



enhancement" 22,23 . We show in the right panel the oscillation probability Pe^, which 
takes relatively large values because of this effect. Here the different curves correspond 
to the matter density profiles in the left panel. Obviously, the two curves are qualitatively 
similar, which means that the three step profile is a reasonable approximation that we apply. 
Note that the parametric resonance cannot be correctly reproduced with a single density 
approximation, leading to qualitatively incorrect conclusions. In the right panel the two 
peaks from the core resonance and mantle resonance can also be seen. From this figure we 



^See discussion in Ref. 29 for a neutrino factory. 
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Figure 2: Comparison between the survival P^^ and appearance P^^ probabilities as a function of 



energy for L = 11810 km and sin^ 26*13 = 0.03, Sep = 0. 

may extract that a detector threshold as low as 2 GeV will be needed to cover the core 
resonance, and a detector threshold as low as 5 GeV to cover the mantle resonance. For a 
threshold much lower than 2 GeV, no qualitatively improvement is expected, whereas for a 
threshold significantly above 10 GeV, no particular sensitivity from the appearance channel 
is expected. 

Having discussed the appearance probability above, let us now focus on the survival proba- 
bility. First, in Fig. [2| we show a comparison between the survival P^^ and appearance Pg^i 
probabilities. It is noteworthy that, for large enough values of sin^2^i3, there are energy 
windows where the appearance probability dominates, and others where the survival proba- 
bility dominates. This effect will be relevant for a neutrino factory if no charge identification 



can be used in the detector 43 and the appearance and survival probabilities cannot be 



distinguished. One can easily see from Fig. ^ that, with sufficient energy resolution, the 



signals from these two probabilities can be separated. Will will quantify in Sec. 5.2 the 
needed energy resolution for this separation. 

For large values of sin^ 2^13, not only the appearance probability contains the information on 
6cp, but also the survival probability P^^ contributes 35 . While a survival probability can. 



in general, not be directly sensitive to the CP violating term oc sin^cp, the dependence on 
cos(5cp is maintained. This may be used to discriminate CP conserving (cos^cp = 1) from 
CP violating (cos^cp ~ 0) solutions. Nevertheless, high statistics are needed to establish 
this signal against the background from the leading (atmospheric) oscillation. We show in 
Fig. [3] the survival probability as a function of energy for different values of sin^ 2^13 (upper 
panel) and different values of 6cp (lower panel). While sin^ 2^13 affects the amplitude of all 
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Figure 3: Disappearance probability (P^^) as a function of energy for different values of sin^ 26*13 (upper 
panel) and different values of ^cp (lower panel). 



oscillation maxima in the same way, a different value of 6cp, compared to the maximally CP 
violating values 7r/2 or 3n/2, leads to a modulation of the oscillation maxima in different 
directions. In principle, sin^ 26'i3 and 6cp can be disentangled by measuring more than one 
oscillation maximum. Again, energy resolution is important for this effect. As indicated 
above, the modulation of the signal is relatively small, meaning that very high statistics are 
required and therefore a large detector may help. 

There is other oscillation physics within the neutrino Standard Model that can be performed 
at very long or core-crossing baselines beyond the sin^ 2^13, 6cp, and the mass hierarchy 
discoveries. For example, such a baseline is helpful (compared to, or combined with a shorter 
baseline) for a measurement of the atmospheric oscillation parameters 44 , deviations from 
maximal atmospheric mixing 44 , determination of the 623 octant 37 , and Sep 45 and 
^13 37 precision measurements. Furthermore, the MSW effect in Earth matter may be 
verified at a high confidence level, even for sin^ 2^13 = 0, by the solar appearance term (the 
fourth term in Eq. (jll)) (46]. Apart from standard oscillation physics, certain classes of new 



physics effects are best tested at such a long baseline, including "non-standard interactions" 
(see Fig. 6 in Ref. [38]). One interesting geophysical application, which we discuss Sec. [7| is 
the determination of the Earth's matter density along the baseline, including that which is 
core-crossing, via the MSW effect 24,37,47 . This apphcation is unique for the setups used 



in this study, since a baseline to the South Pole is the only currently investigated option for 
a core-crossing baseline using a neutrino beam. 
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Label Beam 



Detector 1 Li [km] Detector 2 L2 [km] 



Two-baseline setups for small ^13: 



BBref 


Beta beam (four isotopes) 


WC 


650 


MIND 


7000 


NFref 


Neutrino factory, = 25 GeV 
Beta beam (four isotopes) 


MIND 


7500 


MIND 


4000 


BBbm 


WC 


650 


PINGU 


11810 


NFbm 


Neutrino factory, E^j, = 25 GeV 


MIND 


7500 


PINGU 


11810 


Single- 


■baseline setups for large ^13: 












Beta beam (^^Ne and ''He) 


WC 


650 


n/a 




NF;,f 


Neutrino factory, E^j, = 10 GeV 


MIND 


2200 


n/a 




SB*cf 


LBNE experiment 


WC 


1300 


n/a 




BBU 


Beta beam (^B and ^Li) 


PINGU 


11810 


n/a 




NFbm 


Neutrino factory, E^ = b GeV 


PINGU 


11810 


n/a 




SBbm 


Superbeam (LBNE-like) 


PINGU 


11810 


n/a 





Table 1 : Description and terminology of the reference (subscript "ref " ) and benchmark (subscript "bm" ) 
setups used in this work. 

3 Reference setups 

Here we introduce our reference setups for beta beam, neutrino factory, and superbeam, as 
well as the beam technologies and detector requirements. These reference setups will be 
used later to compare with configurations including the PINGU detector. Note that for 
beta beam and neutrino factory, we define configurations with two baselines (small ^13) and 
one baseline (large ^13). In the following, the two-baseline reference setups will be referred to 
as BBref (beta beam) and NFref (neutrino factory), the single-baseline setups as BB*gf (beta 
beam), NF*^^ (neutrino factory), and SB*j,f (superbeam). The reference setups (subscript 
"ref" ) are summarized in Table [l| 

In all of the following, we perform the numerical simulations using the GLoBES software 



package 48 ,49 . The best-fit oscillation parameters are taken as follows 50 : 

= 33.8° , ^^13 = 5.6° , ^23 = 45° 
Am^^ = 7.7 X lO^^eV^ , Am^^ = 2.4 x 10" W . (3) 

We impose external la errors on Am^i (1%) and 6^12 (1%) and on Amg^ (5%) and 6^23 (5%) 
as conservative estimates for the current measurement uncertainties l50l. We also include a 



5% matter density uncertainty 51 , 52 in the standard performance indicators. Typically, 
we show the results for the simulated normal hierarchy. Note, however, that because of 
the symmetric operation with neutrinos and antineutrinos, the difference for the inverted 
hierarchy is small]^ 

■^The only exception may be the beta beam, where the Ve from '^Li have lower energies than the from 
®B, see below. Here the normal hierarchy performance should be somewhat better. 
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3.1 Beta beam references (BBref, BB*gf) 

Currently discussed configurations of a beta beam produce a neutrino beam by beta decays 
of tlie following proposed isotopes: 

(4) 
(5) 
(6) 
(7) 







I'F + e+ 


+ T^e 


^He 


-> 


^Li + e~ 


+ 




-> 


iBe + 


' + T^e 


lu 




iBe + e" 


' + T^e 



stored in a storage ring with long straight sections (see Fig. 3 in Ref. 53 for such a storage 
ring). The oscillation channels of interest are: 

appearance: z/g — )■ for ^^Ne or *B stored , (8) 

appearance: z/g — )■ for ^'^He or *Li stored . (9) 

Since the beam is flavor-clean, there are no special detector requirements other than {p^) 
flavor identification. This, a beta beam is per se the best candidate for PINGU. 

The isotopes and ^Li have higher endpoint energies than ^^Ne and ^He, which means 
that the same neutrino energies Ei, oc 'jEq can be reached with lower boost factors 7, and 
therefore less acceleration effort. In turn, higher neutrino energies can be reached with the 
same accelerator. Note that the total flux in the detector is proportional to 7^, meaning that 
higher boost factors produce higher event rates. Since beta beams are typically discussed 
in context with the CERN-SPS as accelerator, possibly including upgrades, the maximum 
7 is limited. As a consequence, for small ^13, ^^Ne and ^He are proposed for CP violation 
measurements at shorter baselines, whereas *B and ^Li may be used at a longer baseline for 



determining the mass hierarchy; see Refs. 30 ,54 for the optimization. The chosen detector 



technology for a beta beam typically depends on the neutrino energies: below 1 GeV a 
good separation between the quasi-elastic event from electron and muon neutrinos and an 
efficient background reduction can be reached in a WC detector, a (possibly magnetized) 
iron calorimeter can better identify muon tracks at energies above 1 GeV. 



We use the four-ion setup 30 for reference BBref, which we have reproduced using the 
detector simulation [55,56 . In this setup, ^^Ne and ^He are each accelerated to 7 = 350 and 
detected in a 500 kton WC detector located at Canfranc in Spain with a baseline of 650 km. 
^B and ^Li are accelerated to 7 = 656 and 7 = 390, respectively, and detected at 7000 km 
by a 50 kton magnetized iron detector (MIND), possibly at INO in India. We use a total 
running time of ten years with 10^^ useful ion decays per year in total, effectively 2.5 years 
per each ion species. In practice, the luminosity may be lower, but distributed over a longer 
time, since some of the modes can be run in parallel. This number is ambitious and yet to 



be proven, potentially using the production ring technology 20 . The signal normalization 
errors are 2.5%, with background normalization errors of 5%. The peak neutrino energies 
in this setup are roughly 1.4 GeV for ^^Ne and ^He stored, 9.5 GeV for ^B, and 5.2 GeV for 
^Li. We assume that the straights of the storage ring for the very long baseline are only 60% 
of that of the short baseline storage ring [30j due to the steepness, and this yields 8 ■ 10^^ 
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Figure 4: Discovery reach for CPV, MH and 6*13 at the 3cr confidence level for the reference beta beam 
setups BBrof and BB*^f, where the short and long baseline reaches are shown separately. Note that the 
short baseline of BBicf corresponds to BB*^,^. Here we assume the normal hierarchy. 



useful ion decays in total. For large 6*13, we use a similar setup BB*^,^ with the short baseline 



20 



only, which is close to optimal 31 and also considered within the Euronu design study 
The total number of useful ^B and ^Li decays is the same as for BBj-ef, and no additional 
technology for ^^Ne and ^He is needed, yielding 3 ■ 10^^ useful ion decays in total. 

Shown in Fig. |4] is the performance of the reference beta beam setups as a function of true 
sin^ 2^13 and (5cp (upper row) or fraction of 5cp (lower row) for the CPV, MH, and 6^13 
discovery reaches. In these panels, the parameter of interest will be discovered on the r.h.s. 
of the curve. The thick curves show the two-baseline setup BBj-ef, which have excellent 
discovery reaches for small values of sin^ 26*13. The difference between the one- and two- 
baseline setups can be seen by comparing the dashed and solid curves: the second (longer) 
baseline especially aids degeneracy resolution for small ^13 (left panels) and for the mass 
hierarchy determination for small ^13 (middle panels), whereas it may not be required for 
large 6*13. The large 6*13 reference setup BB*gf (with only one baseline) corresponds to the 
dashed curves. Later we will either replace the long baseline in the two-baseline setup by a 
baseline to PINGU, or the single (short) baseline for large 6*13, and compare the performance 
to these reference setups. 
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3.2 Neutrino factory references (NFref, NF^^f) 

The neutrino factory beam comes from muon decays (muon energy E^) in straight sections 
of a storage ring. Depending on the parent muon, the decay will lead to 

+ : /i" ^ + z/e + z/^ , (10) 
- : /i+ e+ + z/e + z7^ , (11) 

where we have defined the polarities "+" and This implies that for any given muon 
polarity, two different neutrino flavors, one neutrino and one antineutrino, are produced at 
the same time. The oscillation channels of interest are: 

appearance: Ue — )■ i'^ for /i"*" stored (— ) , (12) 

disappearance: -> for /i"*" stored (— ) , (13) 

appearance: Ue for jjT stored (+) , (14) 

disappearance: z/^ — )■ z/^ for /i~ stored (+) . (15) 

Since the different muon polarities can be distinguished by the beam bunching, the neutrinos 
from {u^) appearance and z/^ (z/^) disappearance will be present at the same time in the 
detector, and have to be distinguished. This identification is typically performed with a 
magnetic field that permits the charge identification of the muon produced in a charged 
current interaction. The typical choice for the detector technology is therefore a large 
magnetized iron detector (MIND). Since charge identification is probably unrealistic in 
PINGU, the most practicable other methods to discriminate between the channels may be 



energy resolution 43 ; see discussion in Sec. 2 This, however, is potentially also difficult in 



a detector such as PINGU. 

For a neutrino factory, the small and large 613 setups are specified in the IDS-NF Interim 
Design Report |27j. For small values of ^13, a two-baseline high energy neutrino factory 
(HENF) setup (referred to by NF^ef) with a 100 kt MIND at Li ~ 3000-5000 km and a 50 kt 



72 — 7500 km is proposed, muon energy E^^ = 25 GeV, recently optimized in 29 



38,44 . For large values of 6*13, a single-baseline setup (referred to by NF*^^) with a smaller 



Efj, and shorter baselines is preferable, as proposed as low energy neutrino factory (LENF) 



in Refs. 57 59 . We use the recent data on the MIND detector 17 ,60 ,61 for our simulation. 



This leads to an optimized LENF reference setup with En = 10 GeV and L = 2200 km 29 



for which we use the 100 kt MIND detector. We also include the Ur contamination from the 



decay mode of tau leptons (produced by Ur) into muons in the detector 62,63 , including the 
corresponding migration matrices [29]. Since this background is typically reconstructed at 
lower energies, and may fall into the windows where the appearance information dominates, 
this contamination may be potentially harmful if energy resolution is used for the channel 
discrimination in the detector. We consider the impact of this effect later in the paper. We 
assume a total 10^^ useful muon decays per year, with a total operation time of ten years. 
This number of useful muon decays must be divided among four decay straights (two rings 
with two polarities each) for the HENF, and among two decay straights for the LENF. The 
signal normalization errors are 2.5%, with background normalization errors of 20%. The 
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Figure 5: Discovery reach to CPV, MH, and ^13 at 3a for the reference neutrmo factory setups NFrcf 
and NF*j,f, where the short and long basehne reaches (NFj-ef) are shown separately. Here we assume the 
normal hierarchy. 



peak neutrino energies are about 25 GeV for NFjef and 10 GeV for NF*g£ for muon neutrinos, 
and 17 GeV for NFref and 8 GeV for NF*gf for electron neutrinos. 

Shown in Fig. |5]is the performance of the reference setups. For NFref, the main performance 
for CPV comes from the short baseline, but the long baseline assists degeneracy resolution. 
For the mass hierarchy resolution, the long baseline dominates the performance. For ^13, 
the short baseline performs better, primarily due to the larger detector mass placed there. 
Later, we will replace the long baseline detector by the PINGU detector. For NF*^^, all 
performance indicators can be measured when 613 is large enough, whereas for small ^13 the 
usual degeneracy problems are encountered. In addition, the CPV performance for large 
^13 is better than for NFref. We will compare the PINGU setup for large 613 later to this 
reference setup. 



3.3 Superbeam reference (SB*^^) 

Conventional neutrino beams, such as T2K and MINOS, are produced mainly from the 
decay of pions produced in interactions of a proton beam with a fixed target: 

71+^ + (16) 

vr" ^ fi' + ijf,. (17) 
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Since the muons produced (and also kaons) will further decay, a non-negligible background of 
and is unavoidable. This leads to a systematics limitation for small values of sin^ 2^13, 
which means that superbeams typically have a limited reach in sin^2^i3. It is possible to 



use the off-axis technology 64 to produce a narrow band beam and to over-prop ortionally 
reduce the backgrounds, such as in T2K and NOz/A. In MINOS and LBNE, the beam is 
directly aimed towards the detector, leading to a wide neutrino spectrum. In the LBNE 
case, the spectrum peaks near 2.5 GeV which means that this beam has the ideal energies 
for our purpose. The leading oscillation channels of interest are 

z/e appearance: z/^ — )■ Ve for vr"*" decays, (18) 

disappearance: z/^ — i- for vr''" decays, (19) 

z/g appearance: z/^ — )■ for vr" decays , (20) 

Vfj, disappearance: z/^ — )■ for tt^ decays. (21) 

Thus, for a superbeam detector no particular requirements other than flavor identification 
are necessary, and WC detectors are often proposed. Note that compared to the beta 
beam, the electron flavor must be identified for the appearance channel, which may be 
more challenging than the muon flavor in a detector like PINGU. In either case, since the 
neutrinos carry only a small fraction of the initial proton energy, the intensity of conventional 
beams is rather limited even for high thermal target powers. It therefore remains to be seen 
if such a beam has enough intensity to power a core-crossing baseline. 

We use the 120 GeV configuration of the US-based LBNE project |32j as a prototype. 



referred to as SB*^^. There were two candidates for the far detector: one is a 200 kton 
fiducial mass WC detector; the other is a 34 kton Liquid Argon detector. The baseline is 
1300 km from FNAL to the Homestake mine in South Dakota. We choose LBNE with a 
WC detector as our reference setup for a superbeam experiment, with a total of 120 ■ 10^^ 
(neutrino mode) plus 120 ■ 10^^ (antineutrino mode) protons on target. This corresponds 
to five years of operation in each mode with 2.3 MW target power. For the beam and 



detector simulation, we use the information from the official webpage 32 , and cross checked 



with 65 . The signal and background normalization errors are 5% and 10%, respectively, 
for the disappearance channels, and 1% and 5%, respectively, for the appearance channels. 
We show in Fig. |6] the performance of this reference setup, compared to different detector 
and luminosity configurations. The optimal CPV and MH discovery reaches can only be 
obtained for sin^ 2^^i3 > 0.01 and therefore we only consider this single-baseline setup for 
large 6*13. Later we will compare it to a corresponding PINGU configuration. 



4 Parameterization of the PINGU detector 

Compared to traditional WC detectors, where the photodetectors are placed at the surface 
of a tank, the photodetectors in a conventional neutrino telescope are attached to strings dis- 
tributed within the detection volume. This means that three dimensional space and timing 
information can be used for the event reconstruction, which implies that the event topologies 



in a neutrino telescope are somewhat different from a WC detector 66 . An extraterrestrial 
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Figure 6: A comparison of the discovery reach of CPV, MH and On at Scr for LBNE with different 
himinosities and detectors, including the reference SB*^f. Here we assume the normal hierarchy. 



neutrino detector typically utilizes muon tracks, which are most dominantly produced by 
charged current interactions of muon neutrinos. Other event topologies include electromag- 
netic (from electron neutrinos), hadronic (from tau neutrinos), and neutral current cascades. 
Since it is relatively difficult to distinguish among the different types of cascades, one has to 
accept either large contaminations (flavor-misidentification) for electron and tau neutrinos, 
or apply significant cuts to identify events for which the fiavor can be uniquely assigned. 
In this work, we assume the same fiducial mass for all fiavors as a first approximation. 
This is a reasonable assumption from preliminary fiducial mass computations ^67J, and we 
parameterize the mis-identification rate of the events. In future simulations, the optimal 
balance between cuts and mis-identification has to be determined]^ Other possible event 



topologies include the possibility of a separated tau track at PeV energies 68 . The exact 
computation of event rates (or fiducial masses) for different topologies requires a detailed 
detector simulation of the PINGU detector. 

We use the following properties to parameterize the detector, and we will test the effects of 
each property separately: 

Event mis-identification. We assume that a certain fraction ( "misID" ) of the background 
event rate is mis-identified as signal. For the beta beam, this fraction applies to the neutral 



^Stronger cuts typically imply a lower mis-identification rate coming together with a smaller fiducial 
mass, and therefore lower statistics. A priori it is not clear where the optimum balance between fiducial 
volume (statistics) and mis-identification (systematics) lies in terms of the physics performance, which 
requires some optimization of the cuts. 
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Table 2: Signal and backgrounds used in the parameterization of the PINGU detector. In the GLoBES 
software, the oscillated background rates are calculated, multiplied by "Factor" (and the background sys- 
tematical error parameter), and then added to the signal event rate. Here refers to all flavors, and "CC" 
to charged current interactions. Only one polarity operation is shown for the sake of simplicity. Note that for 
the superbeam, the channels for the other polarity, which is present in the same beam as a contamination, 
are not explicitly shown. 

current and electron neutrino events]^ For the neutrino factory, we take into account that 
the detector is non-magnetized. This means that and events of the same beam polarity 
sum. The mis-identification fraction then applies to neutral current events, z/g, z/g, z/,-, and 
Vr (cascades). In addition, the muon decay mode of the taus ivr and v^) are included 
with a branching fraction of 17% and the corresponding migration matrices, since these 
cannot be disentangled from muon tracks. For the superbeam, we include all neutral and 
charged current backgrounds. However, we distinguish two mis-identification fractions, one 
for muon tracks ("misIDtracks") and one for all the other fiavors ("misID"). Here the logic 
is that muon tracks may be not be as easily confused with other event topologies compared 
to cascades, z.e., misIDtracks < misID. While the disappearance channels are affected by 
misIDtracks, the appearance channels are affected by misID (apart from mis-identified muon 
tracks). The intrinsic beam contaminations are fully taken into account. The signal and 
backgrounds used for the simulation are summarized in Table |2} 

Detection threshold. We impose a sharp threshold on the detection, which corresponds 
to a sharp cutoff of the fiducial mass (energy bins below this cutoff are not taken into 

^The I'e — ^ Tt oscillation channel has the same qualitative dependence on sin^ 20i3 as the v,^ — 
channel, which means that it adds to the sin^ 2^13 signal and hardly affects the measurement. 
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Experiment 


Setups 


misID 


misIDtracks 


Eth 


Energy res. 


Fid. mass 


Beta beam 


BBbm, BB^j^ 


0.001 


n/a 


2 GeV 


50% ■ E 


5 Mt 


Neutrino factory 


NFbm, NF*^ 


0.001 


n/a 


2 GeV 


10% ■ E 


5 Mt 


Superbeam 


SBbm 


0.01 


0.01 


2 GeV 


20% ■ E 


5 Mt 



Table 3: Benchmark values used for the PINGU detector, together with the corresponding setup names. 
The values will be varied and studied individually in the following sections. 

account). From the discussion in Sec. [2} it is clear that such a threshold does not have a 
major effect below 2 GeV. 

Energy resolution We assume an energy resolution AE = t] ■ E, where rj is referred to as 
"energy resolution" . We use a spectral analysis including the energy information, as well as 
we test a total event rate-based analysis. 

Fiducial volume Since the track and shower lengths increase with energy, the fiducial 
volume typically increases with energy. We test different fiducial masses, which we assume 
to be fixed over the whole energy range as a first approximation. We make a comparison to 
preliminary PINGU calculations for the fiducial mass as a function of energy in Appendix [A| 
Note that the main part of this study is independent of a specific detector, meaning it can 
also serve for the optimization of PINGU phase II. 

In the next two sections we vary each of the above parameters to establish the minimal 
requirements for the detector. For each type of experiment, we choose a set of benchmark 
values, listed in Table [3j The logic of the setup names is the same as before: BBbm and 
NFbm correspond to BBj-ef and NFj-ef, respectively. The longer baseline detector is replaced 
by PINGU. BB*^, NF*^, and SB^^^ correspond to BB^^j, NF^^f, and SB^^f, respectively. The 
single baseline is replaced by a baseline to PINGU. The benchmark setups (subscript "bm") 
are summarized in Table [l] Note that some of the values (for example, the energy resolution) 
are not the same for all experiments. These choices reflect that the detector requirements 
for some experiments are more critical than for others. They are chosen to be in the regions 
of highest sensitivity for the detector performance. The values for misID for beta beam and 
neutrino factory imply that muon tracks are identified with a contamination of less than 
0.1% from other flavor^ whereas for the superbeam, such a low mis-identification rate is 
probably unrealistic (at least for electron neutrinos). The benchmark choice for the fiducial 
mass corresponds to the PINGU fiducial mass for z/^ at about 4 GeV from preliminary 
simulations j67j. In either case, we will test the impact of these parameters in the following 
sections. 

The systematic uncertainties are assumed to be the same for the beta beam and neutrino 
factory for consistency: 2.5% signal normalization error and 5% background normalization 
error. For the superbeam, we use 5% and and 10% for signal and background normalization, 
respectively, for the disappearance channels, as for the reference setups. For the appearance 
channels, we adopt a more conservative signal normalization error than for the reference 

®As we will demonstrate later, a mis-identification of 1%, which is perhaps more realistic, will also lead 
to reasonable sensitivities. 
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setup, and we use the same numbers as for the beta beam and neutrino factory. Note 
that these effective systematic uncertainties arise from the beam flux normahzation, the 
near-far extrapolation, and a fiducial mass error, just to name a few examples. Therefore, 
it is a consistent first estimate to use similar numbers for all experiments if the effective 
uncertainties are dominated by the detector. 

5 Requirements for two-baseline setups for small 6is 

Here we discuss the PINGU requirements for two-baseline setups for a beta beam and 
neutrino factory. The second baseline is typically used for the mass hierarchy measurement 
and degeneracy resolution for small values of ^13. Compared to the reference setups in Sec. [3} 
this second baseline is replaced by a core-crossing baseline to PINGU. The performance is 
evaluated compared to the reference setups. 

5.1 Beta beam 

The discovery reach of CPV, MH and ^13 at 3cr for a two-baseline beta beam combination of 
WC detector and PINGU detector is shown in Fig. [7} The solid curves correspond to the 
benchmark values from Table [3j The detector parameters are individually varied in each of 
the rows compared to the benchmark values. The shaded regions show the performance of 
the reference setup, where the right edges of the shadings correspond to the short baseline 
only, and the left edges to the baseline combination. Since the long baseline of the reference 
setup was replaced by the PINGU detector, approaching the right edges of the shadings 
means that the second baseline does not contribute. Exceeding the left edges of the shaded 
regions implies the performance is better than that of the reference setup, which is regarded 
as an optimal setup for each class. A result well within the shaded area means the second 
baseline adds sensitivity, but less than the reference setup. 

Considering each of the detector parameters in turn, the mis-identification fraction (first 
row) should be smaller than about 0.01 to match the original reference setup (left edges 
of shadings), where the MH and ^13 reaches are most sensitive to this parameter. The 
energy threshold (second row) should be at most 5 GeV to match the reference setup, and 
a 2 GeV threshold leads to a significantly better performance than the reference. Beyond 
2 GeV, a further reduction in energy threshold does not show additional improvement. This 
is an expected dependence from the neutrino appearance probability oscillation pattern 
illustrated in Fig. [!} The energy resolution (third row) should be at least about 50% to 
outperform the reference setup. Note, however, that even with a total rate only analysis 
( "rates only" ) the mass hierarchy performance of the reference setup is within reach. Better 
energy resolutions clearly improve the performance, far beyond the reference setup. The 
fiducial mass (lowest row) shows 1 Mt is sufficient to match the performance of the reference 
setup, and larger masses show clear improvement. All things considered, it may be useful 
to improve the mis-identification by appropriate cuts, at the price of a lower fiducial mass 
(efficiency). From Appendix [Al where an energy-dependent fiducial mass is tested, it is 
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Figure 7: Discovery reach for CPV, MH and ^13 at 3cr for a two-baseline beta beam combination of 
WC detector and PINGU detector. The soUd curves correspond to the benchmark values from Table |3] 
i.e., BBbm- The detector parameters are individually varied in each of the rows and compared to the 
benchmark values. The shaded regions show the performance of the reference setup, with the right hand 
edges corresponding to the short baseline scenario, and the left hand edges corresponding to the baseline 
combination. Here we assume the normal hierarchy. 
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Figure 8: Discovery reach for CPV, MH and 6*13 at 3a for different channel configurations of a single- 
baseUne neutrino factory, where a (hypothetical) magnetized detector at the very long PINGU baseline is 
assumed and compared to the short baseline of NFfcf- Values from Table |3] are used. 

clear that the values for threshold and fiducial mass used in Table |3] match such an energy- 
dependent fiducial mass rather well, and are, in this case, possibly conservative. 

In summary, the detector requirements for PINGU as a beta beam detector at a very long 
baseline are modest. The primary issue is a relatively clean fiavor-identification of muon 
neutrinos, which may be achieved by selecting events with some loss of efficiency. On the 
other hand, energy resolution is not so important. 



5.2 Neutrino factory 

A peculiarity for the neutrino factory is the assumption of a non-magnetized detector. This 



means that, for instance, the muon neutrino and antineutrino event rates from Eqs. (12) 



and (13) must be added, and appearance and disappearance channels can only be discrimi- 
nated by energy resolution, as illustrated previously in Fig. |2} However, note that although 
not as efficient, the disappearance channel also has some sensitivity to sin^ 2^13 and ^cp^ 
as discussed in Sec. |2j We show this effect quantitatively in Fig. |8] for a single-basehne 
neutrino factory, where the MIND detector performance at the short baseline is shown for 
comparison. In this figure, we hypothetically assume that the appearance and disappear- 
ance channels could be separated, i.e., the detector would be magnetized. Obviously, the 
very long baseline has an excellent reach in the appearance channels because of the enor- 
mous statistics, and even the disappearance channels can measure all of the performance 
indicators better than the shorter baseline MIND setup. Of course, since the relevant terms 
are sub-leading in the disappearance channels, the performance is significantly worse than 
for the appearance channels. However, if the appearance and disappearance information is 
added, the information on CPV, MH and 613 is significantly reduced compared to either of 
the channels alone (unless ^^13 is very large, where the signal in the disappearance channel 
dominates). This is the situation in a non-magnetized detector, such as PINGU. 

Similar to the beta beam discussion, we show in Fig. [9] the requirements for the long base- 
line non-magnetized PINGU detector as a second baseline detector. As discussed above, 
the most critical requirement for the (non-magnetized) neutrino factory detector is energy 
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Figure 9: Discovery reach for CPV, MH and at 3cr for a two-baseline neutrino factory combination of 
MIND detector and PINGU detector. The soUd curves correspond to the benchmark values from Table |3) 
i.e., NFbm- The detector parameters are individually varied in each of the rows and compared to the 
benchmark values. The shaded regions show the performance of the reference setup, with the right hand 
edges corresponding to the short baseline scenario, and the left hand edges corresponding to the baseline 
combination. Here we assume the normal hierarchy. 
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resolution, where 50%, as used for the beta beam benchmark, is clearly shown to be in- 
sufficient. We therefore apply a smaller value (10%) for the benchmark setup in Table [s] 
On the other hand, at least in principle, the reference setup can be roughly matched for (a 
hypothetical) 5% energy resolution for small 613, whereas already for 20% energy resolution 
a better performance is obtained for CPV for large sin^ 2^13. The dependence on the thresh- 
old (second row) is qualitatively similar to the beta beam, since it depends on the coverage 
of the matter effect peaks. In particular, for small ^13, another critical factor is the fidu- 
cial mass (lower row), where at least 15-20 Mt is needed for a sufficient CPV performance. 
The reason is that for small 613 the appearance signal can only be disentangled from the 
large disappearance background with sufficient statistics. On the other hand, the sin^ 26'i3 
terms in the disappearance channel also be contaminated by a large background from the 
leading oscillation terms. The mis-identification fraction (ffist row) is of minor importance, 
as long as it is, at most, 0.01, because the intrinsic background from the disappearance 
channel is dominant. If an energy-dependent fiducial mass is used, see Appendix |A| our 
estimates in Table |3] are conservative. Compared to the beta beam, the neutrino factory 
spectrum peaks at higher energies, where there actual fiducial mass is larger than 5 Mt. 
Note that we have also included the intrinsic u^- contamination in this analysis. This may 
be potentially harmful since the muon to tau neutrino oscillation sample, which comes from 
a leading oscillation channel driven by the large 623, may be reconstructed at energies where 
the sub-leading appearance channel, driven by the small ^13, dominates. 

In summary we find that the critical requirements for the non-magnetized PINGU as a 
neutrino factory detector are energy resolution and statistics. Due to the higher beam energy 
of a high energy neutrino factory, large fiducial masses may be easily reached. However, we 
note the required energy resolution at a 10% level may be unrealistic for muon tracks in 
such a detector. 



6 Single-baseline setups for large 6is 

In this section we discuss single-baseline setups for large 6*13. We start with the superbeam, 
similar to the previous section, and briefly discuss the optimization of neutrino factory and 
beta beam for large 613. Finally, we compare the precision measurements of the single- 
baseline PINGU setups to the corresponding reference setups. 



6.1 Requirements for the superbeam 



Similar to the previous section, we illustrate the requirements for PINGU as a superbeam 
detector in Fig. 
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As far as the mis-identification is concerned, note that we distinguish 
the identification of muon tracks, as needed for the disappearance channels, and the identi- 
fication of electron neutrinos (cascades), as needed for the appearance channels. Since the 



appearance channels are critical for the measurements in Fig. [10| it can be shown that the 
challenge is the identification of the electron neutrinos. We show in Fig. 10, upper row, the 
corresponding mis-identification fractions. From our tests we find the identification of the 
muon tracks is not critical as long as the corresponding mis-identification fraction misID- 
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Figure 10: Discovery reach for CPV, MH and 6*13 at 3cr for a single-baseline superbeam with core- 
crossing baseline to the PINGU detector. The thick solid curves correspond to the benchmark values from 
Table |3j i.e., SBJ^j^^. The detector parameters are individually varied in the different rows and compared to 
the benchmark values. The thin solid curves show the performance of the reference setup. Here we assume 
the normal hierarchy. 
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tracks < misID, while misIDtracks ^ misID does not significantly improve the results. We 
therefore choose 0.01 for misIDtracks in Fig. 10, showing larger values for misID. Note that 
the discussed mis-identification fractions are substantially larger than the ones for a beta 
beam and neutrino factory, because it could be difficult to assign a specific flavor to the 
cascades. The case for which the flavor of a cascade cannot be assigned corresponds to 



misID=1.0. From Fig. 10 it can be clearly seen that a small mis-identification fraction, 



less than 30%, is necessary to match the performance of the reference beta beam for small 
^13 for CPV. Interestingly, for MH and di^ the performance of the reference superbeam is 
always exceeded even if the cascades cannot be identified. However, the reference CPV 
performance for large 6'i3, i.e., the case for which one would probably build a superbeam, is 
significantly difficult to achieve because of the smaller CPV effects at the long baseline. The 
reason for the "inverted" behavior compared to the reference is that the very long baseline 
is not very sensitive to Sqp and the leading term proportional to sin^ 26'i3 dominates. Note 
that the statistics from the large fiducial mass helps for small ^13. 



As shown in Fig. 10, the dependencies on threshold, energy resolution, and fiducial masses 
are all critical for the CPV measurement. For MH and ^13, typically moderate requirements, 
such as a threshold of 5 GeV and a rates only measurement, are sufficient to match the 
performance of the reference setup. Note, however, that for large sin^ 26'i3 > 0.01, the 
reference superbeam could also measure these parameters for any value of 6cp- In order to 
reach a competitive CPV performance, a relatively large fiducial mass is needed for large 
^^13. However, as seen in Appendix |X| the standard curves already reasonably represent the 
energy-dependent fiducial mass due to the moderately high beam energy. 

In summary, the requirements for PINGU as a single-baseline superbeam detector are 
significant. In particular, for the CPV measurement, good flavor identification of electron 
neutrinos with less than 1% mis-identified events, a large fiducial mass, and a good energy 
resolution are required. Therefore, unless a much larger detector is planned for a second 
phase of PINGU, the experiment may be challenging. The parameter choices in Table |3] 
reflect that result and we note that such a superbeam experiment could serve as a second 
baseline for the mass hierarchy measurement only. 



6.2 Minimal beta beam and neutrino factory setups 

While the superbeam is designed as an instrument for large 613, a neutrino factory and a 
beta beam are typically optimized in a different way. The idea is to propose a "minimal" 
version of these experiments which can perform each of the measurements. In particular, 
only a single-baseline may be used, and lower muon energies (neutrino factory) or 7 (beta 
beam) may be allowed. As with the single-baseline superbeam, the problem is then, of 
course, that the very long baseline is not optimal for CPV measurements. 

We have studied the and 7 dependence for the neutrino factory (NF^^^) and the beta 
beam (BB^j^), respectively. For the beta beam, ®B and ^Li have to be used since the relevant 
oscillation features at the very long baseline otherwise can not be studied fully. We have re- 
scaled the boost factors of both isotopes. For smaller boost factors the CPV discovery reach 
is severely affected, whereas larger boost factors may not be possible even at an upgraded 
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Figure 11: Fraction of ^cp for which CPV will we observed (3(t) as a function of i?^ for sin^ 2^13 = 10~^. 
Here the energy resolution is varied, and the effect of the Vt backgrounds is shown. See Table [3] for the rest 
of the benchmark values used. The label NF*gf refers to the low energy (i?^ = 10 GeV) reference version of 
a neutrino factory with MIND at 2200 km. Here we assume the normal hierarchy. 



CERN-SPS. Therefore, we use the earher proposed 7 = 656 and 7 = 390 for and ^Li, 
respectively. 

For a neutrino factory, the mass hierarchy and Q\z can typically be measured at the very 



long baseline for large 6'i3 and the critical factor is the CPV discovery reach. In Fig. 11 the 
dependence of the fraction of 5cp for the CPV discovery as a function of E^^ for a large value 
of sin^ 2^13 = 0.1 is shown, which is close to the current bound. For the reference setup 
(MIND@2200 km) it can be clearly seen that around 10 GeV is optimal, for which CPV 
can be measured for at least 65% of all (5cp. It is instructive to study the curve "PINGU 
without r" with 10% energy resolution and without the intrinsic r contamination (z/^ or 
producing tau leptons, which decay into muons). This curve monotonously increases with 
Ey^ and plateaus between about 7 and 18 GeV. This behavior can be understood in terms 
of Fig. [2] the Pg/i-dominated range near 7.5 GeV is covered by the spectrum, where the 
discovery reach will be saturated. It is interesting to note the behavior if the r contamination 
is included. The curve is similar to the one without r contamination until about 6 GeV, 
since the r production threshold is not significantly exceeded, and the performance then 
decreases. This is due to 17% of the from the leading z/^ — )■ channel being mis- 
identified as v^J^ adding to the Ve — ?■ appearance channel in the non-magnetized detector. 
Note that these events are reconstructed at lower energies, which means that the ranges in 
Fig. [2| where the appearance information dominates, are affected. Only if Ey^ approaches 
25 GeV, the performance does increase again. Therefore, for the LENF to PINGU, we 
choose Ey = 5 GeV, where the r contamination issue can be completely avoided and the 
comparison with the reference setup shows a similar performance. As it can be seen from 
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Figure 12: The precision measurements of CP phase i^cp and sin^ 29i3 for three single-basehne neutrino 
experiments: Beta Beam (BB), Neutrino Factory (NF), and SuperBeam (SB). The contours represent the 
la, 2a and 3a confidence levels (2 d.o.f.). Filled contours represent the PINGU benchmark setups and 
unfilled contours the reference setups. The crosses mark the best fit value of sin^ 26*13 and 6cp- Here 
we assume the normal (true) hierarchy. The inverted (fit) hierarchy solution can be ruled out by the 
experiments. 

the curves with different energy resolutions, this choice is at the lower end of the allowed 
range, but closer to the original idea of the LENF. From the 15% energy resolution curve it 
can be also seen that the r contamination must be avoided in that case; here the power of 
the energy resolution to discriminate different channels is severely limited by the different 
reconstructed energies of the events from r decays. 



6.3 Precision measurements with single-baseline setups 



After the recent release of long-baseline electron neutrino appearance searches at the T2K 
and MINOS experiments, a global analysis hints at a la level to the Dirac CP violating 
phase pJO]. We also consider the best fit values of sin^ ^13 = 0.013 and 6cp = — O.GItt for 
the normal hierarchy and study precision measurements at three neutrino experiments with 
different single-baseline beams in Fig. 
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The single-baseline reference setups are defined 
in Sec. |3| whereas the PINGU versions are described earlier in this section. Note that there 
is only a fit for one hierarchy in the figure, because the inverted hierarchy solution can be 
excluded. 



From Fig. [T2| it can be seen that the PINGU setup can not only compete with neutrino 
factory and superbeam reference, but, in fact, has performance which is somewhat better. 
For the neutrino factory, this could be see already in Fig. [9} where the long baseline helped 
for Sep if sin^ 2^13 was large. For the superbeam, the PINGU setup in Fig. 12 (right panel) 



appears to be better, which is opposite the behavior seen in Fig. 10, where, for large ^13, the 
reference setup is better. The reason for this is the different dependence on 6cp- The CPV 
performance roughly corresponds to the precision at 6cp = or vr, but the precision shown 

the performance of 



in Fig. 


12 


Fig. 


10 


is 



-O.Glvr. We confirm that for Scp 



Sqp- For a relatively narrow beam spectrum in particular, depending on the baseline, the 
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Figure 13: Relative precision of the Earth's core density measurement as a function of the true value of 
sin^ 2^13 for three neutrino experiments with different beams at the Icr, 2a and Scr confidence levels. Here 
we assume the normal hierarchy with Jcp = 0. 



sin^cp- or cos (5cp-term can be emphasized, see Eq. ([T]), leading to a strong dependence on 

For tlie beta beam, it was already obvious from Fig. [7] that the CPV does not significantly 
benefit from the very long baseline. This explains why the CPV is significantly worse than 
for the shorter baseline reference in Fig. 12 Compared to neutrino factory and superbeam, 
the main issue here is the relatively asymmetric operation with and Ve with very different 
peak energies, whereas the muon energies for i/g and z/e at the neutrino factory are exactly 
alike. This is particularly true since the Vf. spectrum peaks at relatively low energies and 
the disappearance channel is missing. 

In summary, for large 6'i3, PINGU may be an interesting option to receive a beam sent 
from the LENF or a seperbeam. We note, however, that the values required for the detector 
performance are relatively challenging, see Table [3j For instance, a very good energy resolu- 
tion is required for both setups. On the other hand, the beta beam, for which the detector 
performance is of secondary importance, suffers from constraints on the boost factors from 
the CERN-SPS. A balanced green-field scenario {i.e., a setup with a new accelerator) may, 
however, have a better performance. 



7 Measurement of the Earth's core density 

A feature of the PINGU setups discussed in this work is that they cross the Earth's outer 



core. It may be possible to use this to determine the core density by the MSW effect 24 
Here we discuss the quantitative performance for such a measurement of the two-baseline 
(beta beam, neutrino factory) or single-basehne (superbeam) benchmark setups, discussed 
earlier in this work. The density profile of the Earth for the long baseline is approximated 
by the three-step profile in the left panel of Fig. [Tj For the sake of simplicity, we assume that 
the Earth's mantle density is relatively well known from seismic wave geophysics, and that 



25 





la 


TV iT 1 C 1 1 1 ' 1 

Measurement oi the core density e 
2a 


3(T 


BBbm 


+U.U/D3 

-0.0461 


+0.053/ 
-0.1030 


+0.0807 
-0.1380 


NFbn, 


+0.0058 
-0.0059 


+0.0117 
-0.0120 


+0.0174 
-0.0181 




+0.0046 
-0.0045 


+0.0091 
-0.0089 


+0.0137 
-0.0132 



Table 4: The relative precision of core density measurement e = [p^^l — p^^f-'^) / p^^f^"^ at different confidence 
levels assuming sin^ 20i3 = 10~^. 



the core density, i.e. the higher density in Fig. [T| is to be measured]^ In this measurement, 
the oscillation parameters are marginalized over, which means that the uncertainties on the 
oscillation parameters, as obtained from the experiment itself, are included. 



We show in Fig. 13 the relative matter density precision for the (outer) core density mea- 
surement as a function of the (true) sin^2^i3, assuming that the chemical composition is 
known. For large sin^2^i3, the information in the appearance channel is typically clean 
enough for a relatively precise measurement (the backgrounds from the second to fourth 
terms in Eq. ([l]) are small). For intermediate values of sin^ 2^13, correlations and degenera- 
cies prohibit the measurement of the density. For small values of sin^ 26'i3, the measurement 
is dominated by the solar term; c/., last term in Eq. ([T|. 

For the neutrino factory, the core density can be determined to an accuracy of up to 0.5% 
at la for sin^ 2^13 > 10~^. For small values of sin^ 26'i3, however, the solar term cannot be 
easily extracted because of the disappearance channel contamination in the non-magnetized 
detector. For the superbeam, the performance is similar to that the neutrino factory, and 



to predictions from past studies 37,47 . Precise numbers at different confidence levels are 
found in Table |4| assuming that sin 2^13 = 10~^. As a major difference of superbeam 
and beta beam compared to the neutrino factory, the matter density can be also measured 
relatively well for small values of sin^ 2^13 because the appearance and disappearance chan- 
nels can be separated. However, for large sin^26'i3, the statistics of the beta beam are not 
comparable to that of neutrino factory and superbeam, for the reasons mentioned in the 
previous subsection. 



8 Summary and conclusions 

A future upgrade of the IceCube/DeepCore detector at the South Pole, known as PINGU 
("Precision IceCube Next Generation Upgrade"), would have a lower threshold via denser 
photomultiplier spacing, making it, in principle, an interesting candidate for a neutrino 
beam. While the baseline from any of the major accelerator laboratories in the Northern 
hemisphere is not optimal for CP violation studies, the huge increase of the detector fiducial 
mass, compared to conventional detectors, may compensate. We have tested the require- 



^See Ref. 



37 for the test of more complicated hypotheses. 



26 



ments for such a detector, which is currently in the design phase. As a starting point, we 
have used beta beam, neutrino factory, and superbeam setups which are accepted to be 
optimized for either small (undiscovered) or large (discovered) sin^2^i3. For the small 613 
setups, we have replaced the second "magic" baseline by a baseline to PINGU. For large 
6*13, confirming the recent results from T2K and Double Chooz, we have discussed "mini- 
mal" single-baseline versions. In all cases, we have compared the performance quantitatively 
to the known optimized reference setups in order to determine the minimal requirements 
for the PINGU detector. As main performance indicators, we have used the conventional 
discovery reaches for CP violation, mass hierarchy, and ^13. For large 613, we have also 
shown the precision for the ^13 and 6cp measurements. 

From the point of view of neutrino oscillations, the baseline to PINGU is interesting since 
it crosses the Earth's (outer) core. This implies the appearance channel is parametrically 
enhanced, see Fig. [l| with resonance energies corresponding to the Earth's mantle and 
core densities. This has immediate consequences for the detection threshold of the PINGU 
detector and the required beam energies: the core resonance requires a measurement between 
2 and 5 GeV, and the mantle resonance between 5 and 10 GeV. Thresholds lower than 
2 GeV or significantly above 10 GeV minimally impact the results. The chosen beams in 
this study take this energy range into account, and it is coincidentally a range where the 
PINGU fiducial mass is large. 

For small values of ^13, a scenario where the recent T2K and Double Chooz results are not 
confirmed and PINGU has replaced the magic baseline detector, we have identified that 
the flavor-clean beta beam imposes the least requirements to the detector. Most important 
is a relatively clean flavor-identification of muon neutrinos. Using muon tracks, this may 
be achieved by applying strict cuts at the expense of fiducial mass. Energy resolution 
is of secondary importance. The neutrino factory, on the other hand, requires channel 
identification by energy resolution if the detector is non-magnetized. The required energy 
resolution of about 10% E may prove challenging for short muon tracks in the PINGU 
detector. 

For large ^13 (^13 discovered), PINGU has been studied as a single-baseline standalone 
detector of a long-baseline experiment. For a superbeam, using the LBNE beam as a 
guide, the requirements are high for CPV measurements: flavor identification of electron 
neutrinos with less than 1% mis-identified events from other flavors, a large fiducial mass, 
and a quality energy resolution are needed. Unless a much larger detector is planned for 
a phase II of PINGU or the detector is used for the mass hierarchy measurement only, 
the experiment may be challenging. For the beta beam, the CPV has found to be strongly 
constrained by the limited boost factors in an upgraded CERN-SPS. However, a green-field 
scenario not restricted to this particular accelerator may have better performance. For the 
neutrino factory, energy resolution is a major challenge, similar to the small sin^ 2^13 case. 
However, small muon energies as low as E^j^ 5 GeV are preferable, where the intrinsic r 
contamination (muon tracks from tau decays) can be avoided. This problem has turned out 
to be especially severe in a non-magnetized detector, where the secondary muons (from the 
tau decays into muons) from the leading channel — )■ u^- may be reconstructed in energy 
windows where the appearance channel is to be isolated by energy resolution. Such a low 
muon energy may be especially interesting in connection with recent discussions on a low 



27 



luminosity and low energy initial stage neutrino factory. 

We have also studied the potential measurement of the Earth's outer core density using the 
MSW effect, which is special for PINGU setups. We have found that, depending on setup 
and sin^ 26'i3, a precision on the outer core density as high as 0.5% at 3a can be reached for 
a superbeam and a neutrino factory. 

In conclusion, for small 6'i3, the beta beam may be the most promising beam option for 
the PINGU detector, since it is a flavor-clean beam with moderate requirements for the 
detector. For large 613, a superbeam or neutrino factory might be an interesting alternative, 
provided that an excellent energy resolution can be obtained and other detector requirements 
can be met. In particular, a low energy neutrino factory may benefit from the parametric 
enhancement and the core resonance energy of a few GeV. In addition, the intrinsic tau 
contamination can be avoided by the low beam energy. Finally, note that we have only 
studied parameter sets inspired by phase I of such a proposal, whereas it is an option to 
increase the fiducial mass in a later phase. In this case, one should further optimize for 
fiducial mass in the range between 2 and 10 GeV instead of a lower threshold from the 
point of view of long-baseline physics. 
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A Energy-dependent effective fiducial volume 

The effective area or fiducial volume of a neutrino telescope depends on the neutrino energy, 
since a larger energy produces a longer muon track or cascade, and thus deposits more 
light. In the extreme case of low energies, only the volume close to the photomultipliers will 
contribute to the fiducial volume, whereas for very high energies even light from outside of 
the actual detector volume may be detected. In the main text of this study, we parameterize 
the fiducial volume by a fixed value with a threshold. Here we show the effects of an energy- 



dependent fiducial volume generated for a preliminary simulation of PINGU phase I 67 
Note that while the conclusions in the main text do not depend on the actual geometry of 
the detector, the conclusions in this Appendix depend on the fiducial volumes. 



We show the fiducial volume as a function of energy for the different flavors in Fig. 14 This 
volume is defined as 

Events triggered ^1 

Veff = 7; X Generation volume , (22) 

Events generated 
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Figure 14: The effective fiducial volume of a preliminary PINGU configuration as a function of the 
incoming neutrino energy for three flavor neutrinos. The dashed-dotted line is the constant fiducial volume 
chosen for the benchmark setups in the main text. The peak energies of the different neutrino sources are 
marked. 




sin^29i3 sin'26i3 sin^26i3 



Figure 15: A comparison of the discovery reach of CPV, MH and 613 at 3(7 for the benchmark setups 
(c/., Table [3]) with or without the effective volume Kff of Fig. 14 The benchmark setups without energy- 
dependent Vas are represented by the solid curves, while those with energy-dependent Vcs are represented 
by the dashed curves. Here we assume the normal hierarchy. 



i.e., it has to be interpreted as a function of ttie reconstructed neutrino energy. Since 
it includes charged current and neutral current interactions, and it does not include cuts 
required to improve the flavor identification, and can only be regarded as an estimate. We 



also show in Fig. 14 the peak energies of the different instruments {E^ = 25 GeV for the 
neutrino factory), and the value for the fiducial volume used as standard value in the main 
text (5 Mt). One can easily extract from this figure (middle panel) that our estimate is a 
good match for the low energy neutrino factory denoted by NF* and beta beam denoted 
by BB(^B/'^Li). However, the fiducial volume for the high energy neutrino factory, denoted 
by NF, is actually significantly underestimated and the effective volume for the superbeam 
(left panel) is slightly overestimated. Note that the beam has a high energy tail, which may 
be useful in parts of the parameter space. From the right panel (z^t), one can also extract 
that our z/^ contaminations, which are based on the same fiducial volume as z/g and in 
the main text, are somewhat overestimated. 
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Finally, in Fig. 15 , a quantitative comparison between the setups based on fiducial volume 



and threshold in Table |3] and those based on Fig. 14 are shown in (two-baseline setups for 
beta beam and neutrino factory used). This figure confirms the previous estimates for the 
full simulation. 
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